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Abstract 

The relations among the occupation number of the lowest natural orbital (ONLNO), momentum distributions (MD) 
and off-diagonal long-range element (ODLRE) of the reduced single-particle density matrix (RSPDM) are studied 
while Tonks-Girardeau gas in one dimensional periodic potential is in the ground state. For A^-body systems of 
large enough, RSPDM and its lowest natural orbital do not vary with N in overlapped areas in commensurate and 
incommensurate cases correspondingly. In commensurate case, the ODLRE is exponential attenuation with N, which 
results in that the ONLNO and MD are invariant with N. While in contrast, in incommensurate case, the off-diagonal 
elements are inversely proportional to y/N, which results in the different behavior of the ONLNO and MD. 
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1. Introduction 

With the development of realizing quasi-one dimen- 
sional cold atomic systems in experiment||jj,|2t], one di- 
mensional systems which could be exactly solved have 
again drawn people's attention. One of these is Tonks- 
Girardeau Jl, 01 (TG) gas which consists of bosons with 
"impenetrable" core repulsive interactions. Actually, at 
low temperatures and densities, Bose system will act as 
a TG system|5l|3,01. TG gas is first a Bose gas and to 
some extent it also displays some properties just like a 
Fermi gasfli]. 

As the TG gas model can be exactly solved via Fermi - 
Bose mapping[4 ll0lll 111 which does not depend on the 
external potential, many models with different exter- 
nal potentials have been well studied theoretically. 5- 
split harmonic potential model was studied by J. Goold 
and Th. Busch lfl2ll . They showed that the scaling of 
the occupation number depends on whether one has an 
even or odd number of particles. Kronig-Penney poten- 
tial model was investigated by Lin, etc. il 311 and Wei, 
etc. |3, who pointed that there are two different phases 
in periodic potentials. One is Mott insulator phase that 
the ratio of the number of bosons to wells N/M is an 
integer, the other is boson conductor phase which N/M 
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is a fractional number. The two phases displayed en- 
tirely different physical properties II 1 411 . Cosine poten- 
tial model was considered by Z. L. Wang and A. M. 
Wang[9]. There are two advantages of choosing cosine 
potential. First, it is easier to prepare ground state in 
Mott insulator phase. Second, it is easier to prepare 
cosine potential than 6 potential and it is convenient to 
control amplitude B and frequency at. In this work, we 
continue studying cosine potential. 

Many properties of the TG gas and its correspond- 
ing Fermi gas are always identical!! H, [Till , such as 
the average value of particle coordinates, potential en- 
ergies, system total momenta, single-particle densities 
and pair distribution functions. But other properties 
are quite different, such as RSPDM and MD functions 

Before R. Pezer 

and H. Buljan||8|] introduced a method for calculating 
the RSPDM of a TG gas, it is difficult to calculate the 
RSPDM and the co rresp onding MD for large number 
svstemsfl5, "l6L 17, 18, Ts|. In this manuscript, we 
mainly use Pezer' s method[8] to study non-local prop- 
erties of the TG gas. 

Generally speaking, there are three discrimination 
methods in judging whether a boson system occurs 
Bose-Einstein Condensate (BEC), such as the exhibi- 
tion of off-diagonal long-range order (ODLRO), the 
macroscopic occupation numbers of the lowest natural 
orbital and the zero-momentum state|20ll2lLl22fl. These 
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physical quantities are all derived from RSPDM, so they 
are not independent and we believe that they have some 
relations among them. Through the studies of proper- 
ties of TG gas in periodic potential, we have found the 
relations and deduced one form in the text. 

This work is organized as follows. In Sec.2, we intro- 
duce the model Hamiltonian and its computation pro- 
cesses. In Sec. 3, we study the ground state properties 
of the TG gas by numerical calculation and theoretical 
explanation. In Sec .4, brief conclusions are given. 

2. Model Hamiltonian and wave functions 

In a TG system, the boson is assumed to have an 
"impenetrable" hard core characterized by a radius of 
a. From Girardeau's work, with the hard core radius 
a — » 0, the interparticle interaction is given by 



U(x,,x.) 




(1) 



Such an interparticle interaction could be represented 
by the following subsidiary condition on the wave func- 
tion iff : 



i/fix, 



, x N , = if x. = x r l^i<j^N. (2) 



With Girardeau's Bose-Fermi mapping [4], the wave 
function of TG system is given by wave function of 
noninteracting spinless Fermions multiply an antisym- 
metric factor. 



in which 



,x s ,t)=A(x 1 



• x N )i/^(x. 



,x N ,t), (3) 



A(X U - ■ ■ ,X N ) = Y\ SgnO; - Xj), 



(4) 



i>] 
1, 



-1, 



x > 0, 
x < 0. 



Now we study the case that magnetized bosons in 
an external magnetic field B(x) = -Bcos(2wx). With- 
out considering the interparticle interactions, the total 
Hamiltonian is written as 



"=2 



2m dx 2 



in which 



V(x) = jiB cos(2<jjx), 



L L 

2 <X< 2 



(5) 



(6) 



We suppose L — Mn/aj, where M is an integer, /j. 
is the magnetic moment of atom and to is the circular 
frequency. For the sake of simplicity, we assume that 
both N and M are odd (As a result of the introduction of 
A(Xj , • • • , x N ), the wavefunction ifr of TG system is pe- 
riodic if is odd while otherwise antiperiodicfl I23I0 . 
The single-particle Schrodinger equation inie [—k , |] 
is written as 



h 2 d 2 

-— — +fxB cos(2a>x) 
Lm ox L 



(p m (x) = E a (p m (x). (7) 



„ , . muB 
Substitute z = wi, q = — — - and A 



tion Q, then it becomes 

d 2 <p(z) 



h 2 u> 2 



2mE 

—r— into equa- 



dz 2 



+ [A- 2q cos(2z)]ifim(z) = 0, 



(8) 



Z G 



ojL ojL 

Using Bloch's theorem and periodic boundary condi- 
tions (p m (z) = (p m (z + &*L), fmiz) is periodic with period 
n, so we can expand it in Fourier series: 



<P> 



= exp(i'vz) ^ c n exp(;2nz), 



(9) 
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v = — , /e0,±l,- 
M 1 



M-l 



Then substitute (|9]l into ([8]), we obtain an infinite sym- 
metric tridiagonal matrix equation which is the eigen 
equation. 
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= A 



(10) 



By truncating the matrix in each direction (centered 
at the smallest diagonal element v 2 ) at sufficiently large 



2 



dimensions, approximations to the desired eigenvalues 
and eigenfunctions can be obtained to any desired pre- 
cision. For more detailed computation process, please 
refer to the references!!, B [jj B . 



The wave functions of the TG system are given by the 
Slater determinant 

iff = A(* lt • ■ ■ ,xJ-i=DeC . =1 [p m ( Xj )] . (11) 



Nl 

When the temperature is zero, the system will be in 
the ground state — N particles in the N lowest eigen- 
states, respectively. In this work, we mainly consider 
the ground state, so we should choose the lowest N 
single-particle eigenstates in Eq. (fTTl i. 

3. Relations Among Some Properties 

It is easy to prove that single particle density and pah- 
distribution functions are always the same for TG gas 
and its mapping spinless Fermi gasfl,S[l3, 23]. Here 
we just consider properties that differ from the two sys- 
tems, such as RSPDM, MD, natural orbital and its oc- 
cupation numbers. 
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Figure 1 : (Color online) RSPDM p(.v, x' ) of the TG gas in external 
periodic magnetic field for different particles in commensurate and 
incommensurate cases, q = 0.1, 01 = 1. The abscissa axis x and 
the ordinate axis x' are in units of n. (a) N = 15, N/M = 1; (b) 
N = 15, N/M = 3/5; (c) N = 51, N/M = 1; (d) N = 51,N/M = 3/5. 



3.1. Reduced Single-particle Density Matrix 

The reduced single-particle density matrix, with nor- 
malization J p(x, x) dx = N, is given by 



p(x, X 



')=N J f{x, 



•> x N )if/(x', x 2 ,- ■ ■, x N )dx 2 - ■ -dx N . 

(12) 

The RSPDM, which expresses self-correlation, re- 
flects the probability that having found a particle at po- 
sition x and at the same time detecting the particle at 
position x' . Many important observable physical quan- 
tities are defined by the RSPDM, such as the occupation 
numbers of natural orbital and the MD, but it is hard to 
calculate p(x, x') by the definition of equation (12), even 
if numerical calculate. Fortunately, it can be expressed 
in terms of the dynamically evolving single-particle ba- 
sis raised by R. Pezer and H. Buljan|8]: 

N 

p(x, x) = ^ ip*(x)A u (x, x')ipj(x'). (13) 

The N xN matrix A(x, x') = {Ajj(x, x')} is 

A(x,x') = (p-yDetP, (14) 

where the entries of the matrix P are Pjj(x, x') = Sy - 
2 f * dy tp*(y)ipj(y), and we have assumed x < x' without 
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Figure 2: (Color online) Off-diagonal elements p M (0, x) in commen- 
surate case, N/M = 1. q = 0.1, 0) = 1. From the higher to the lower: 
N= 15,25,35,57, 123. 



loss of generality. In addition, the RSPDM in periodic 
potential satisfies IU4I1 



p(x + T, x + T) — p(x, x'), 



(15) 



where T is the period of external periodic potential, and 
in this work T — n/u). With this formula, we can calcu- 
late considerable large systems of the TG gas. 

Figure Q] displays contour plot of RSPDM for both 
the Mott insulator and boson conductor phase. For both 
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Figure 3: (Color online) Off-diagonal elements p c (0, x) in incommen- 
surate case, N/M = 3/5. q = 0.1, 01 = 1. From the higher to the 
lower: N = 21, 33, 51, 75, 123. 



Figure 4: (Color online) Off-diagonal elements - In p M (0, L/2) as a 
function of N in commensurate case, N = M, q = OA, at = 1. 



the two phases in large systems, p(x, x') does not vary 
with N in overlapped areas for fixed N/M. In order to 
illustrate this view, we have plotted p M (0, x) for N - M 
in figure [2] and p c (0, x) for N/M = 3/5 in figure [3] 
(The reason of choosing p(0, x) is to avoid the bound- 
ary effect[Ili], which results from the periodic bound- 
ary conditions), where the subscripts M and C stand for 
Mott insulator phase and boson conductor phase respec- 
tively. As the system is in Mott-insulator phase, all par- 
ticles tend to be localized and have little correlations in 
long distance. Therefore, p(0, x) decreases quickly in 
Mott-insulator phase while relatively smooth in boson 
conductor phase. 

3.2. Off-diagonal Elements 

ODLRO was suggested by Yang [22] as the type of 
ordering for superfluids, such as BEC and electron pairs 
in superconductors. In the thermodynamic limit and \x- 
x'\ — » oo, if the off-diagonal element p(x, x') = 0, then 
ODLRO is not present and there is no BEC; Otherwise, 
ODLRO is present and there is BECl|20i[2ll[22|. 

Considering the periodic boundary conditions, we as- 
sume p(0, L/2) to be a typical value of off-diagonal el- 
ements in the thermodynamic limit. The off-diagonal 
elements p M (Q, L/2) and p c (0, L/2) as a function of par- 
ticle number N are shown in figure [4] and ^respectively. 
Based on these numerical data, empirical correlation on 
ODLRE was obtained via regression. The results are as 
follows, 

p M (0,L/2) = C 1 exp(-/fiA0, (16a) 

p c (0,L/2) = C 2 N- l/2 , (16b) 

where C\, C2 and /ci are constants and they vary with q. 
From equation ( TTBT l. we know the off-diagonal elements 
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Figure 5 : (Color online) Off-diagonal elements - In p c (0, L/2) as a 
function of InN in incommensurate case, N/M = 3/5, q = 0.1, oj = 
1. 



p(0, L/2) are zero in both the commensurate and incom- 
mensurate cases in the thermodynamic limit. Therefore, 
there are no ODLRO and BEC in both commensurate 
and incommensurate cases of TG gas in the external pe- 
riodic potential. 

3.3. Occupation Numbers and Natural Orbital 

The occupation numbers and the natural orbital are 
defined as 



dx'p(x, x')(/)i(x') — Ai<pi(x), 



(17) 



where Aj represents the occupation number of the natu- 
ral orbital <pj and A, = N. For simplicity, we label the 
eigenvalues A; in a descending order: Aq > A\ > A2 > 
■ ■ ■ . The corresponding state 0o of is the condensate 
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state of bosons. For the ground state, <po is the lowest 
natural orbital and the corresponding occupation num- 
ber Aq is the largest. For large systems, from equation 
(fTol l we know the off-diagonal elements with long dis- 
tances have few contributions to 0o> so the scaled low- 
est natural orbital VA^>o(jc) are the same in overlapped 
areas for different N in commensurate and incommen- 
surate cases separately. The lowest natural orbital is pe- 
riodic with period T = n/co, which is the same as the 
external periodic potential. We have plotted two peri- 
ods in figure [6] 
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Figure 6: (Color online) Scaled lowest natural orbital <po(x) for N = 57 
and N = 123 in commensurate case N = M and incommensurate case 
N/M = 3/5, q = 0.1, OJ = 1. 



nothing to do with x, so the result of the second term in 
equation (TT~8b should not include x. We set x — with- 
out loss of generality. Using the mean value theorem for 
integrals, the second term becomes 



<h(N, 



— f 

V,0) Jn 



< 1 + ff)2 



0-*>3 



dx'p(0, x')<p (N, x') * Cp (0,L/2). 

(19) 

where C is a constant. Rewriting equation ( TT~8T >. we ob- 
tain the relation between the occupation number of the 
natural orbital and the off-diagonal elements, 



AAo 
AN 



A (N)-A (N-l) = Cp(0,L/2). 



(20) 



In the thermodynamic limit and regard N as continuous, 
equation ( fZOl i becomes 



dAo 
dN 



= Cp(0, CO), 



(21a) 



or 



~dN 



C lim lim therm p(x, x), (21b) 



where "lim therm" means "thermodynamic limit". 
From equation ( f2TT > we know that the behavior of the 
largest eigenvalue Aq and the off-diagonal elements 
p(x, x') as \x - x'\ — > co are one-to-one correspondence. 



In a macroscopic system, the presence or absence 
of BEC is determined by the behavior of p(x, x') as 
\x - x'\ — > co or the largest eigenvalue Aq of p(x, x') 
proportional lEolEllE^ to N. Both the occupation num- 
ber Ao and the off-diagonal element p(0, L/2) are deter- 
mined by p(x, x'), hence they must have some deep rela- 
tions. For large systems, it is convenient to use (po(N, x) 
and Aq(N) to represent the lowest natural orbital and its 
occupation number of A^-body system. Then, 



If",, 

: 1 /at \ dx P( x > x )^o(A 7 , x ) 

1 f(!-w)T 
=— — dx'p(x,x')MN,x') (18) 

MN,x) J- ( i-i)i 

—— dxp(x,x )4> (N,x ). 



From the definition of /I,-, we know the first term of the 
right hand in equation ( fT"8b is just Aq(N - 1), where we 
have used the invariance of p(x, x') and the accompany- 
ing V^V0o(A f , x) for fixed N/M. The value of A (N) has 
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Figure 7: (Color online) ONLNO as a function of N in commensurate 
case, q = 0. 1, ai = 1. 

Substitute equation ( fTol l into ( 121 al l, we can get the ex- 
pressions of Aqm an d ^oc m commensurate and incom- 
mensurate cases separately, 



^om = a[l - expf-pW)] + y, 



(22a) 
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Figure 8: (Color online) ONLNO as a function of N in incommensu- 
rate case N/M = 3/5, q = 0.1, to = 1. 
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Figure 9: (Color online) Normalized MD, q = 0.1, to = 1. (a) and 
(b): Normalized MD for different N in commensurate and incommen- 
surate (N/M = 3/5) cases, (c) and (d): n(0) as a function of N in 
commensurate and incommensurate (N/M = 3/5) cases. 



A oc = a' y/N + ■ 



(22b) 



where a, /3, y, a' and y' are constants, y and y' are the 
correction terms resulted from the periodic boundary 
conditions. In both the two phases, Aqm or Aoc equals to 
1 while N - 1, so we can choose y = 1 -a[l -exp(-yS)] 
and y' = 1 - a'. We have plotted Aqm and Aqc and their 
fitted curves with equation (1221 in figure Q and [8] They 
clearly show that the numerical result is consistent with 
our theoretical equations. 

Generally speaking, a! > 1 in one dimensional pe- 
riodic TG systems. In this case, we can prove that 
equation (I22bb is equivalent to Girardeau's result lfl6ll 
Aoc ~ where cr is a constant that slightly larger 
than 0.5 and it gradually diminishes to 0.5 in the ther- 
modynamic limit ! 1411 . 

3.4. Momentum Distribution 

The normalized momentum distribution, related to 
the RSPDM, is defined as 



n(k) 



= \ p{x,x')e- ik(x - x,) dxdx', (23) 

2nN J 



with the normalization 



n(k) dk = 1 . 



Actually, MD is the Fourier transformation of the 
RSPDM and they are one-to-one correspondence. 
Therefore, the behavior of MD will be determined by 
the behavior of the elements of RSPDM. We have plot- 
ted n(k) in commensurate and incommensurate cases in 
figure [9] (a) and (b) separately, in which the platforms 



will diminish as /V increases and they will vanish in the 
thermodynamic limit. 

The quantity n(k) is the probability density of par- 
ticles, on the average, in the momentum k state. If a 
momentum state ko has a macroscopic occupation such 
that n(ko) °c N, then the system will exhibit BECli27tl. 
In this work, it is clear to see that ko = is the state 
that has the largest occupation from figure|9](a) and (b). 
With the method used in equation ([T8l >. we can prove 
that n(0) has the similar relationships as equation (fJTJ, 



dn(Q) 
dN 



= C'p(0,oo), 



(24a) 



dn(Q) 

— — = C lim lim therm p(x, x ), (24b) 

dN \JC-X?\-HX> 

where C is an constant. 

Substitute equation ( fT6b into ( 124al ). n M (0) and « c (0) 
in the commensurate and incommensurate cases can be 
expressed as 

n M (0) = ^[l-exp(-7 7 /V)]+r, (25a) 

n c (0) = £' Vat + t', (25b) 

where 77, t and t' are constants. r\ equals to f3 

while the external periodic potential and particles are 
the same, where /3 is a parameter in equation d22l) . t 
and t' are the correction terms resulted from the pe- 
riodic boundary conditions. n M (0), n c (Q) and their fit- 
ted curves with dZSb are shown in figure |9](c) and (d). 
It is easily to see that the number of particles occupy- 
ing zero-momentum state n(0) is not proportional to N, 
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Figure 10: (Color online) Aq as a function of n(0), in commensurate 
(a) and incommensurate (b) cases, q = 0.1, ai = 1. 



which is another evidence that there are no BEC in TG 
gas in periodic potential. 

The relations between the ONLNO Ao and «(0) can 
be obtained by divide equation (12 lab and ( I24ab . 



dn(0) 



= C, 



so 



Cn(0), 



(26) 



(27) 



where we have used the fact that Ao = n(0) = as 
N = 0. The relations between Ao and n(0) are plotted in 
figure[10](a) and (b) in commensurate and incommensu- 
rate case separately. Since the momentum spectrum can 
be calculated by just Fourier transforming the natural 
orbital, it seems to me that the linear relation between 
Aq and n(0) is not surprising. 

4. Discussions and Conclusions 

From equations ( 121 al l and d24ab . we know that Aq and 
n(0) are proportional to N if p(0, oo) is a constant and 
nonzero, so BEC occurs in this case. Otherwise Aq and 
n(0) are not proportional to N and there is no BEC. T his 
inference is accord with the former workKC 



We have studied the ONLNO, MD and off-diagonal 
elements of RSPDM in the ground state of TG gas in 
an external periodic magnetic field. The three quanti- 
ties are depicted as characteristics of judging whether 
a Bose system exhibits BEC or not. We have found 
the relationships among them, which are expressed as 
equations d2Tl i, d24b and d26l ). In the special case of one 
dimensional periodic TG system and in the thermody- 
namic limit, the ONLNO Aq and MD n(0) are propor- 
tional to A/ in commensurate case while to y/N in in- 
commensurate case. 
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